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Theoretical predictions for the BCS-BEC crossover of trapped Fermi atoms are compared with re- 
cent experimental results for the density profiles of 6 Li. The calculations rest on a single theoretical 
approach that includes pairing fluctuations beyond mean field. Excellent agreement with experi- 
mental results is obtained. Theoretical predictions for the zero-temperature chemical potential and 
gap at the unitarity limit are also found to compare extremely well with Quantum Monte Carlo 
simulations and with recent experimental results. 
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The original theoretical ideas behind the crossover, 
from the Bardeen-Cooper-Schrieffer (BCS) superconduc- 
tivity with largely overlapping Cooper pairs to the Bose- 
Einstein condensation (BEC) of dilute bosons, date back 
to the pioneering work by Eagles for low-carrier doped 
superconductors 1 . Later seminal papers by Leggett and 
by Nozieres and Schmitt-Rink have provided a general 
framework for the BCS-BEC crossover, both at zero tem- 
perature in the superfluid phase 2 and at finite tempera- 
ture in the normal phase 3 . The experimental motivations 
to these studies came from the condensation of excitons 
in solids , pairing in nuclei 5 , and pseudogap in high- 
temperature superconductors 6 . No direct quantitative 
comparison between theory and experiments has, how- 
ever, been possible so far, owing essentially to the large 
number of degrees of freedoms present in these systems. 

Recent experimental advances on the condensation 
of ultracold trapped Fermi atoms 7 make it now possi- 
ble to compare theoretical predictions for the BCS-BEC 
crossover with experimental results. In these systems, 
the use of a tunable Fano-Feshbach (FF) resonance 8 pro- 
vides the fermionic attraction that triggers pairing, thus 
enabling one to pass with continuity from weak (BCS) to 
strong (BEC) coupling across the crossover region. In ad- 
dition, for broad enough resonances, the scattering length 
a,F appears to be the only relevant quantity entering the 
many-body Hamiltonian of the interacting Fermi atoms 9 . 
Experiments on ultracold Fermi atoms thus constitute an 
ideal testing ground for theories which describe the pro- 
gressive quenching of the fermionic degrees of freedom 
into composite bosons. 

Comparison of theory with experiments is more inter- 
esting for the BCS-BEC crossover with Fermi atoms than 
for the BEC with Bose atoms. This is because the di- 
luteness condition for Bose gases makes it appropriate to 
describe the bosonic condensate by the Gross-Pitaevskii 
equation 10 and the excitations above it by the Bogoli- 
ubov approximation 11 . For the BCS-BEC crossover, on 
the other hand, many-body approximations can be con- 
trolled only on the two (BCS and BEC) sides of the 
crossover, where the diluteness condition holds for the 
gas of fermions and composite bosons, respectively. A 



small parameter to control the many-body approxima- 
tions is instead lacking in the intermediate (crossover) 
region where the scattering length <if diverges. 

Any sensible theoretical approach to crossover phe- 
nomena sets up a single theory that recovers controlled 
approximations on both sides of the crossover and pro- 
vides a continuous evolution between them. It is then 
of particular relevance that for the BCS-BEC crossover 
the results for the intermediate-coupling regime can be 
compared with accurate results from Quantum Monte 
Carlo (QMC) simulations 12 performed at the unitarity 
limit where aj? diverges. This provides a further strin- 
gent test on the validity of a BCS-BEC crossover theory. 

Characteristic of any BCS-BEC crossover theory is to 
provide two coupled equations for the order parameter A 
and the chemical potential \x (the latter being strongly 
renormalized from one limit to the other) . In the present 
theory for the trapped gas, these equations are obtained 
by considering a local-density approximation to the the- 
ory of Ref. 13. The overall chemical potential is replaced 
whenever it occurs by the local quantity /x(r) = fx — V(r) 
that includes the trapping potential V(r) at position r, 
as discussed in Ref. 14. In Ref. 13, the theory devel- 
oped by Popov for a weakly-interacting (dilute) super- 
fluid Fermi gas 15 was extended as to include the effects 
of the collective Bogoliubov- Anderson mode 16 in the di- 
agonal fermionic self-energy. This was obtained by con- 
sidering the self-energy (with Nambu notation) 

Sn(fc) = -| j^T^MSiiiq-k) (1) 
E 12 (fc) = -A, (2) 

where r u (q) = Xu(-q)/bcu(<l)Xn(-Q) -Xuiq) 2 } is the 
normal pair propagator with 
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In these expressions, q = (q, Q v ) and k = (k, uj n ) (q and 
k are wave vectors, ilj, and u) n bosonic and fermionic 
Matsubara frequencies, respectively), m is the fermion 
mass, T the temperature, and Qij (i,j = 1,2) are BCS 
Green's functions. (We set S=fcs = l throughout.) 

In the BCS limit, the quantity kpap (ap < 0) iden- 
tifies the small parameter needed to control the approx- 
imations. Here, the Fermi wave vector Uf is obtained 
by setting k F /(2m) equal to the Fermi energy Ep of the 
non-interacting Fermi gas. The scattering length is 
varied in experiments with trapped Fermi atoms from 
negative to positive values across the FF resonance where 
ap diverges. It was predicted 14,17 that the crossover 
from fermionic to bosonic behavior in a trap occurs, in 
practice, within the narrow range —1 < (/cf^f) 1 < 1. 
By our generalization of the Popov fermionic approxi- 
mation, in the strong-coupling limit of the fermionic at- 
traction the same fermionic theory is able to describe 
a dilute system of composite bosons within the Bogoli- 
ubov approximation 11 . In this limit, the small parameter 
kpap (ap > 0) corresponds to the "gas parameter" of the 
dilute Bose gas. 

Experiments with trapped Fermi atoms are usually 
made using anisotropic harmonic trapping potentials, 
having different frequencies w x , uj y , and lo z for the 
three axis. In this case, the Fermi energy Ep equals 
C&NijJxUJyUJz) 1 ! 3 where TV is the total number of Fermi 
atoms in the trap. Within a local-density approxima- 
tion, the anisotropic problem can be readily mapped onto 
a corresponding isotropic problem with the same value 
of the Fermi energy. Detailed experimental axial den- 
sity profiles for 6 Li Fermi atoms were reported in Ref. 18 
across the intermediate (crossover) region of interest. In 
Ref. 18, uj x = uj y — w r and lu z = Xuj r with A <C 1. The 
experimental values of u> r , u> z , ap, and TV determine the 
value of (kpap)~ l , to be used for comparison with the 
theoretical calculations. In addition, for the experiments 
reported in Ref. 18 the temperature is estimated to be 
much smaller than Ep. Correspondingly, the theoretical 
calculations can be performed at T = 0. 

Figure 1 shows the comparison between our theoretical 
predictions for the axial density profiles in the crossover 
region and the data reported in Fig. 4 of Ref. 18, for three 
different values of the magnetic field B tuning the FF 
resonance. The intermediate value B = 850G (about) 
corresponds to the position of the resonance at which 
aF diverges. For the two other values B — 882G and 
B = 809G, the values of the scattering length are esti- 
mated to be ap = —1.8 x 10 4 ao and of = 8.5 x 10 3 ao, re- 
spectively, where ao is the Bohr radius. These values are 
extracted from the data reported in Fig. 3(a) of Ref. 18. 
The value of kp is obtained with uj r /(2ir) = 640Hz and 
uj z /(2it) = ^/0.6 x B[G] + 32Hz, and depends weakly on 
the estimated value of TV which is affected by the largest 
experimental uncertainty (of the order of 50%). 
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FIG. 1. Comparison between experimental and theoretical 
axial density profiles. Experimental data from Ref. 18 (dots) 
are shown for three diffeent values of the magnetic field B tun- 
ing the FF resonance. Theoretical results at T = obtained 
by our theory (full lines) and by BCS mean field (dashed lines) 
are shown for the corresponding couplings (kpap)" 1 given in 
the text. The upper (lower) panel refers to the estimated 
number of atoms TV = 4 x 10 5 (TV = 2.3 x 10 5 ). 

This uncertainty affects only the absolute scale of the 
experimental density profiles but not their shape. The 
upper panel of Fig. 1 corresponds to the estimated value 
TV = 4 x 10 5 given in Ref. 18. With the above values 
of the frequencies, the coupling (kpap) 1 is completely 
determined to be —0.20, 0.00, and 0.43 from left to right, 
in the order. Our corresponding theoretical predictions 
(full lines) compare well with the experimental results 
(dots) in all three cases. The agreement between theory 
and experiment becomes almost perfect when considering 
the smaller value TV = 2.3 x 10 5 (a value which is within 
the bounds of the experimental uncertainty) , as shown in 
the lower panel of Fig. 1. In this case, (kpap) -1 equals 
-0.22, 0.00, and 0.47 from left to right, in the order 19 . 

Figure 1 reports also the results of the BCS mean-field 
calculation (dashed lines). It is evident that the agree- 
ment between theory and experiment is improved by our 
theory (full lines) which includes quantum pairing fluc- 
tuations beyond mean field. Nevertheless, the results of 
the BCS mean field appear to be reasonably good in com- 
parison with experiment. This comparison thus verifies a 
long-standing theoretical expectation 2 that at zero tem- 
perature the BCS mean field should constitute a reason- 
able approximation for the whole BCS-BEC crossover. 

Experimental data over a wider range of (/cf^f) -1 are 
also available, corresponding to an extended variation 
of the magnetic field about the FF resonance. Figure 
2 shows the comparison between our T = theoreti- 
cal predictions for the normalized root-mean square axial 
radius (triangles) and the data reported in Fig. 3(c) of 
Ref. 18 (dots), versus the coupling (kpap)" 1 in the range 
— 1.5 < (fcj?ai?) _1 < 1.5 spanning the whole crossover. 
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FIG. 2. Comparison between experimental (dots) and the- 
oretical (triangles) normalized root-mean square axial ra- 
dius across the crossover regime. Experimental data are 
taken from Ref. 18. The values of (^fcif) -1 and of the 
non-interacting root-mean square axial radius used also for 
the experimental data are obtained with N = 2.3 x 10 5 . 

The root-mean square axial radius has been normal- 
ized to its expression y/ Ep / '(4mu^) for non-interacting 
fermions. The values of (kpciF^ 1 and of the non- 
interacting root-mean square axial radius (needed to 
compare the experimental data with our theoretical re- 
sults) have been obtained with N — 2.3 x 10 5 , corre- 
sponding to the lower panel of Fig. 1. The agreement 
between theory and experiment is remarkably accurate 
even over this wide coupling range. 

As mentioned already, comparison between a crossover 
theory and the experimental data is most com- 
pelling in the intermediate-coupling regime, due to the 
lack of a small parameter to control the many-body 
approximations 22 when the scattering length ap diverges. 
Special features further occur with a local theory at the 
unitarity limit (kpap)^ 1 = 0. In particular, this limit 
corresponds in the isotropic case to the universal density 
profile (with r =| r |) 
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which depends on the single parameter (3. This parame- 
ter was introduced experimentally in Ref. 23, and can be 
extracted theoretically from the ratio n/Ep = 1 + (3 as 
obtained at the unitarity limit for the homogeneous sys- 
tem. For given values of the harmonic frequencies and of 
(3, the density profile in Eq. (5) thus depends only on the 
total number of atoms N via Ep. 

Recent T = QMC simulations 12 have yielded the 
value (3 = —0.56 ± 0.01. Our theory gives (3 = —0.545, in 
excellent agreement with these simulations. By contrast, 
T = BCS mean field gives (3 = -0.41. Our theoreti- 
cal value of j3 is also fully consistent with recent experi- 
mental data, which yield the values — 0.68^° io (Ref- 18) 
and —0.55 ± 0.10 (Ref. 20). The agreement between 



our theory and QMC simulations is further confirmed 
by comparing the values of the order parameter A at 
the unitarity limit for the homogeneous case. Our the- 
ory gives A/Ep — 0.53, while the QMC simulations of 
Ref. 12 provide the estimated value A/Ep = 0.54. It is 
then apparent that our crossover theory, which captures 
the essential physics on the two sides of the crossover, is 
also able to provide quantitatively accurate results in the 
intermediate-coupling regime, where no small parameter 
can be identified to control the many-body problem. 

It is worth commenting that the number of atoms 
N = 2.3 x 10 5 , utilized by our calculations in the lower 
panel of Fig. 1 as well as in Fig. 2, was determined by 
equating the experimental value for the root-mean square 
axial radius at the unitarity limit given in Ref. 18, to the 
corresponding value obtained from the universal profile 
of Eq. (5) with (3 taken from the QMC simulations of 
Ref. 12. We have already seen from Fig. 1 (lower panel) 
and Fig. 2 that this value of N makes the agreement 
between our theory and the experimental data remark- 
ably good. This value of iV also improves considerably 
the agreement between the experimental value (130nK) 
as given in Ref. 18 for the bosonic chemical potential 
fiB deep in the bosonic regime, and the corresponding 
theoretical estimate (125nK) obtained from the Gross- 
Pitaevskii theory 10 with the value as = 0.6af for the 
dimer-dimer scattering length as ■ If one would instead 
take N = 4.0 x 10 5 (as in the upper panel of Fig. 1), a 
larger value (155n.fr) would result for the theoretical es- 
timate of fls- 

In the present approach, the system of fermions maps 
in strong coupling onto a system of composite bosons 
(dimers), which are described by the Bogoliubov the- 
ory with a dimer-dimer scattering length a^/ciF = 2. 
This value corresponds to the Born approximation for the 
dimer-dimer scattering 25 . It has been shown 25 that in- 
clusion of higher-order dimer-dimer scattering processes 
within the many-body diagrammatic theory decreases 
this value to as/ap = 0.75. A scattering approach to 
the four-fermion problem 24 has given instead the value 
clb/o^f = 0.6 quoted above. The difference between 
the values 0.75 and 0.6 originates from additional dia- 
grams which were not considered in Ref. 25 albeit they 
survive in the zero-density limit. It is interesting that 
the difference between the values 2 and 0.6 — 0.75 can 
be appreciated, by analyzing the experimental results 
for the low-temperature axial density profile reported in 
Fig. 1(b) of Ref. 18 deep in the BEC region (correspond- 
ing to (kpap) 1 = 5.44). Figure 3 compares the experi- 
mental results (dots) directly with the predictions of the 
Bogoliubov theory, obtained with the alternative values 
cib/clf = 2 (dash-dotted line), as/aF = 0.75 (dashed 
line), and as/ap — 0.6 (full line). It is evident from this 
figure that a correct treatment of the dimer-dimer scat- 
tering improves the comparison with experimental data 
in this extreme BEC limit, where the dimer-dimer scat- 
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tering length is the only relevant interaction parameter. 
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FIG. 3. Comparison between experimental results for the 
axial density profile (dots) and the predictions of Bogoli- 
ubov theory, obtained with os/of = 2 (dash-dotted line), 
Ob/of = 0.75 (dashed line), and asjaF = 0.6 (full line). Ex- 
perimental data were obtained in Ref. 18 deep in the BEC re- 
gion, corresponding to (/cfOf) 1 = 5.44 when N = 2.3 x 10 5 . 

The analysis of the many-body diagrams made in Ref. 25 
suggests, however, that inclusion of diagrams beyond 
the Born approximation for the dimer-dimer scatter- 
ing should become immaterial when approaching the 
crossover region, as these diagrams correspond to high- 
order pairing-fluctuation processes of the Fermi system. 

We comment finally that in our approach we have 
not introduced independent bosonic entities besides the 
fermions. The fermions themselves bind, in fact, into 
composite bosons even at finite temperature and density, 
provided their mutual attraction is sufficiently strong. 
Correspondingly, the many-body problem contains the 
two-fermion bound state via an effective single-channel 
problem with scattering length op, which can be identi- 
fied directly with the scattering length varied experimen- 
tally via the FF resonance. The excellent agreement with 
the experimental data shown above demonstrates that 
this single-channel treatment is appropriate to describe 
the BCS-BEC crossover for trapped 6 Li Fermi atoms. 

We are indebted to R. Grimm, J.E. Thomas, and K. 
Levin for discussions, and to R. Grimm for providing us 
with the data files of the figures in Ref. 18. This work was 
partially supported by the Italian MIUR (contract Cofin- 
2003 "Complex Systems and Many-Body Problems"). 
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